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ABSTRACT 

The growth rate of solar activity in the early phase of a solar cycle has been known to be well 
correlated with the subsequent amplitude (solar maximum). It provides very useful information 
for a new solar cycle as its variation reflects the temporal evolution of the dynamic process of 
solar magnetic activities from the initial phase to the peak phase of the cycle. The correlation 
coefficient between the solar maximum (/3 a ) and the rising rate (/3 a ) at Am months after the solar 
minimum (i? m i n ) is studied and shown to increase as the cycle progresses with an inflection point 
(r = 0.83) at about Am — 20 months. The prediction error of -R max based on /3 a is found within 
estimation at the 90% level of confidence and the relative prediction error will be less than 20% 
when Am > 20. From the above relationship, the current cycle (24) is preliminarily predicted to 
peak around October, 2013 with a size of i? m ax = 84 ± 33 at the 90% level of confidence. 

Subject headings: solar physics; solar activity; sun spots; solar cycles 



1. Introduction 

The Waldmeicr effect that stronger cycles tend 
to rise faster is a well known fact in solar activity 
(Waldmeierlll939HHathaway et al.ll2002i: |Pu et al 



2()()9ali. The growth rate of solar activity (i? z ) in 
the early phase of a solar cycle, defined as the ra- 
tio of a given increment (Ai? z = 20) between two 
certain levels (R z i = 30 and i? z2 = 50) over the 
corresponding elapsed time (At), was found to be 
highly corr elated (r > 0.8) with the s ubsequent 
amplitude ( Cameron fc Schusslerl [20081 . There- 
fore, the strength of a new cycle should be ratio- 
nally predicted by the above relation. The prob- 
lem is that if and at what month after the start of 
a new cycle the strength of the cycle can be well 
estimated from the early information of the cycle. 

This paper studies the variation of the correla- 
tion between the maximum amplitude (R max ) of 
a solar cycle and the rising rate (/3 a ) as a function 
of Am months entering the cycle and analyzes the 
predictive power of j3 a on i? max in order to find 
out at what month i? max can be well estimated 
by /3 a . The results are shown in the following sec- 
tion. j3 a is defined as the ratio of the increment 



of R z from the minimum (i? m i n ) over the elapsed 
time (Am months). The temporal variation in the 
correlation coefficient (r) between i? max and /3 a is 
analyzed in Scction [2.1i showing that r is very low 
near the initial phase (r < 0.5 if Am < 10) and 
significant only at a few months after the start of 
the cycle (r > 0.8 if Am > 19). The predictive 
power of j3 a on i? max as the cycle progresses is an- 
alyzed in Section 12.21 indicating that the relative 
prediction error of i? max is very small for almost all 
Am in some cycles and smaller than 20% at some 
(about twenty) months after the start in other cy- 
cles. The peak size and its timing of cycle 24 are 
estimated in Section [31 followed by conclusions in 
Section [H 

2. Data and Analysis 

The data used in the present study are the 
smoothed monthly mean international sunspot 
number (i? z )0 from July, 1749 to February, 2011. 
The rising rate is defined as the ratio, /3 a = 
(R z (Am) — i? m i n )/Am, of the increment of R z 



1 http: / /www. ngdc.noaa.gov/stp/spaccwcather. html 
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Table 1: Parameters and results 



n 


R • 


R 


T a /? a (27,n) 


E (E t ) 


cc (cc t ) 


1 


8.4 


86.5 


76 


0.74 








2 


11.2 


115.8 


40 


2.23 








3 


7.2 


158.5 


35 


3.94 








4 


9.5 


141.2 


41 


3.64 








5 


3.2 


49.2 


82 


0.45 








6 


0.0 


48.7 


70 


0.20 








7 


0.1 


71.5 


79 


0.63 


0.18(0.80) 


— 0.15( 


—0.98) 


8 


7.3 


146.9 


40 


3.58 


0.15(0.27) 


_0.90( 


—0.98) 


9 


10.6 


131.9 


55 


1.24 


0.32(0.31) 


+0.39( 


—0.88) 


10 


3.2 


98.0 


50 


1.54 


0.13(0.47) 


+0.11( 


—0.98) 


11 


5.2 


140.3 


41 


2.49 


0.09(0.32) 


— 0.32( 


—0.98) 


12 


2.2 


74.6 


60 


1.76 


0.39(0.57) 


-0.35( 


-0.98) 


13 


5.0 


87.9 


47 


2.34 


0.30(0.49) 


+0.18( 


-0.98) 


14 


2.7 


64.2 


49 


1.42 


0.43(0.67) 


-0.65( 


-0.96) 


15 


1.5 


105.4 


49 


1.94 


0.09(0.41) 


-0.01( 


-0.93) 


16 


5.6 


78.1 


57 


1.93 


0.22(0.52) 


-0.30( 


-0.95) 


17 


3.5 


119.2 


43 


1.99 


0.15(0.35) 


-0.29( 


-0.96) 


18 


7.7 


151.8 


39 


2.73 


0.21(0.28) 


-0.81( 


-0.98) 


19 


3.4 


201.3 


47 


5.26 


0.19(0.20) 


-0.96( 


-0.95) 


20 


9.6 


110.6 


49 


2.42 


0.06(0.38) 


-0.22( 


-0.97) 


21 


12.2 


164.5 


44 


3.14 


0.25(0.26) 


-0.93( 


-0.98) 


22 


12.3 


158.5 


34 


4.64 


0.10(0.31) 


+0.65( 


-0.98) 


23 


8.0 


120.8 


47 


2.21 


0.09(0.35) 


-0.82( 


-0.98) 


24 


1.7 


?84.0 ?59 


1.17 








av. 


5.9 119.1 


49 


2.43 


0.20(0.41) 


-0.31( 


-0.96) 



from the minimum (i? m ; n ) over the elapsed time 
(Am) from the start of the cycle. The rising rate 
is computed for each cycle n and each Am, de- 
noted by /3 a (Am,n). The parameters are listed 
in Table [1] in which i? m j n and -R max are the mini- 
mum and maximum amplitudes of the solar cycle, 
respectively; T a is the rise time from minimum to 
maximum; /3 a (27, n) is the value of /3 a (Am,n) at 
the current state Am = 27; other parameters will 
be described later; and the last row indicates the 
relevant averages over cycles n = 7-23. 

2.1. The variation in the correlation be- 
tween i? max and /3 a as the cycle pro- 
gresses 

Figure [T] illustrates the variation in the cor- 
relation coefficient (r) between i? max (n) and 
j3 a (Am,n) for the cycles in which Tjji) > Am 
at a given Am (using only the data in the rising 
phases). One can see that r varies with the pro- 
gression of the cycle (Am). A steady increasing 
trend is shown in r since Am = 6: r increases 
from about 0.33 at Am = 6 to about 0.83 at 
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Fig. 1. — Correlation coefficient (r) between i? max 
and /3 a as a function of Am. 

Am = 20 at a high speed, and increases at a 
smaller speed since then, showing an inflection 
point at Am = 20 months, r(20) = 0.83. Near 
the initial phase (Am < 10) the correlation is not 
strong (r < 0.5). The correlation coefficient be- 
tween -R max and /3 a is high enough at Am = 19 
(r > 0.81) months entering the solar cycle. At the 
current state (Am = 27), r(27) = 0.88. 

2.2. The predictive power of /3 a on i? max 

In order to test the predictive power of /3 a on 
i? mm at different Am, we use only the data up to 
cycle (n — 1) to predict i? max for cycle n. For 
a given Am, we calculate the linear regression 
equation of i? max (i) against /3 a (Am,i) for cycles 
i = 1, 2, • • • , n — 1 in the form of 

flmax = A + Bp a , (1) 

and the standard deviation <t(Ato, n — 1) used 
to estimate the uncertainty of the prediction 
of i? max (^)- Then, substituting the value of 
/3 a (Am, n) into this equation, the i? max value 
for cycle n can be predicted, which is denoted 
by R p (Am,n). Figure [2] shows the results for 
the recent nine cycles n = 15-23: R p (Am,n) 
(black solid line) together with error bars t r (n — 
l)tr(Am,n — 1); i? max (n) (black horizontal long- 
dashed line), the actual relative prediction error 
(red dotted), 

£(Am,n)= |flp(Am ' n) ~ E "- (n)l ; (2) 
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the estimated relative prediction error (green 
dashed), 

„ , . . t r (n — llfffAm, n — 1) 

Et Am, n = ^ -fA> '- , 3 

-R max (n) 

where t T (n — 1) is the t- value at the 90% level of 
confidence in a student's t-distribution for ri[ — 
(n — 3) degrees of freedom; and the correlation co- 
efficient between R max (i) and /3 a (Am, i) (r, blue 
dash-dotted line, multiplied by 100 to be indicated 
by the right hand scale). The numbers in the fig- 
ure (cc) denote the correlation coefficients between 
E (E t ) and r. 

Figure [^a) illustrates the results for cycle 15. 
It is seen that i? ma x (black horizontal long-dashed 
line) is all within the error bars of R p at the 90% 
level of confidence (vertical lines), E < E t , and 
E < 20% (red horizontal long-dashed line) when 
Am > 3, although there are some fluctuations in 
both E and E t . The anti-correlation coefficient be- 
tween E t and r is very strong, cc t = —0.93, imply- 
ing that the higher the correlation coefficient (r) 
between i? max and /3 a , the smaller the estimated 
relative prediction error (Et) from the extrapo- 
lation of the relationship between i? max and /3 a . 
However, this is only an estimate in theory rather 
than in practice. In fact, the correlation coefficient 
between E and r is almost zero, cc = —0.01, im- 
plying that the actual relative prediction error is 
almost uncorrelated to r. That is to say it is un- 
certain whether a more (less) accurate prediction 
corresponds to a higher (lower) correlation. 

In Figure EJd) for cycle 18, we test only the re- 
sults for the rising phase Am < T a (= 39). One 
can see that both E and E t decrease as Am in- 
creases. The anti-correlation coefficient between 
Et and r is also very strong, cc t = —0.98. The 
correlation coefficient between E and r is highly 
negative, cc = —0.81, implying that a more (less) 
accurate prediction can be obtained from a higher 
(lower) correlation in this case. In cycle 18, i? max 
is always within the error bars of R p at the 90% 
level of confidence, E < E t . In addition, E < 20% 
when Am > 26. 

Figure [^g) shows the results for cycle 21: 
Am < T a (= 44). The results are similar to 
those in Fig. [2jd): both E and E t decrease as 
Am increases; the anti-correlation coefficient be- 
tween E (E t ) and r is very strong, cc = —0.93 
(cc t = -0.98). At a small Am, E is large: E > E t 



if Am < 17. As the cycle progresses, E becomes 
smaller: E < E t when Am > 18; E < 20% when 
Am > 26. 

The results in other cycles are similar to those 
above. The main conclusions in Fig. [2] may be 
summarized as follows. 

1- Rmax (black horizontal long-dashed) is usu- 
ally within the error bars of R p (black solid) 
at the 90% level of confidence (vertical lines) , 
apart from cycles 19 and 21 when Am < 16 
and Am < 17, respectively; 

2. the estimated relative prediction error E t 
(green dashed) tends to decrease as Am in- 
creases; E t is highly anti-correlated with the 
correlation coefficient r (blue dash-dotted), 
cc t » -0.96; 

3. the actual relative prediction error E (red 
dotted) varies with some fluctuations and 
tends to decrease as Am increases in some 
cycles (for n = 16-19, 21, 23); 

4. there is no established relationship between 
E and r, such as \cc\ > 0.8 in cycles n = 18, 
19, 21 and 23, while cc| < 0.3 in cycles n = 
15, 16, 17 and 20, and cc is even positive in 
cycle 22 (0.65); 

5. E < E t for all Am in cycles n = 15-18, 20, 
22 and 23; E < E t for Am > 17 in cycle 19, 
and for Am > 18 in cycle 21; 

6. E < 20% (red horizontal long-dashed) since 
a very few months entering the cycle (Am > 
3) in cycles 15, 20 and 22; E < 20% for 
Am > 34 in cycle 16, for Am > 21 in cycle 
17, for Am > 26 in cycle 18, for Am > 17 
in cycle 19, for Am > 26 in cycle 21; and for 
Am > 10 in cycle 23, 

In summary, /3 a behaved very well in predicting 
the subsequent -R max : (i) the actual prediction er- 
ror (known only when the cycle is over) is usually 
within estimation since about twenty months en- 
tering the cycle, \R p (Am, n) — R max (n)\ < t T (n — 
l)cr(Am, Ti — 1); (ii) the relative prediction error is 
usually less than 20% since about twenty months 
into the cycle; and (iii) E tends to decrease as 
the cycle progresses. In some cycles (n — 15, 20 
and 22, see Figures EJa), (f) and (h)), the rela- 
tive prediction error is very small (E < 10%) at 
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Fig. 2. — Predictions (i? p , black solid line, left hand scale) of i? max (black horizontal long-dashed line) 
together with error bars (t T (n — l)a(Am, n — 1), vertical line) for cycles 15 - 23 (panels (a) - (i), respectively), 
the actual relative prediction error (E, red dotted, right hand scale), the estimated value (E t , green dashed), 
and the correlation coefficient between i? max and /3 a (r, blue dash-dotted). The numbers in the figure (cc) 
denote the correlation coefficients between E (E t ) and r. 



a small Am even if the correlation coefficient is 
low (r < 0.5). Similar conclusions can also be ob- 
tained in other cycles (not shown): E < E t at 
about Am > 20; E t is highly anti-correlated with 
r (cc t ); while there is no established relationship 
between E and r (cc). The results of E (E t ) and 
cc (cct) in cycles 7-23 are shown in Table [TJ and 
the relevant averages are indicated by the last row: 
< E >= 0.20 (< E t >= 0.41) and < cc >= -0.31 
(< cct >= —0.96), where E is the average over 
Am in a solar cycle and < E > represents the av- 
erage over cycles 7-23. Therefore, a higher (lower) 
correlation coefficient does not n ecessarily yield a 
more (less) accurate predictio n (|Du et al . 2009b; 
Du fc Wandl2011atlDull2011ah . 



3. Prediction i? max for Cycle 24 

Now, we employ the above technique to pre- 
dict the peak size of cycle n — 24. The results 
are shown in Fig. [3j R p (solid) is the predicted 
i? max (24) and r (dotted) is the correlation coef- 
ficient between i? max (i) and /3 a (Am, i) for cycles 
i = 1,2, ••■,23 at a given Am. It is seen that 
R p does not vary significantly with Am. At the 
current state (Am = 27), the correlation coeffi- 
cient between i? max and (3 a is r(27) = 0.88, and 
the regression equation of i? max against /3 a is 



Rn 



52.1 + 27.2/3 a , 



(4) 



with a standard deviation of a = 19.2. Sub- 
stituting the current value of /3 a (27, 24) = 1.17 
into this equation, the peak sunspot number for 
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r(27)=0.88 



R(24)=84.0 
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Prediction of i? max (24) as a function of 



the ongoing cycle (24) is predicted as i? p (24) = 
84±i r (23)cr = (84±33) (asterisk), where t r (23) = 
1.721 is the t-value at the 90% level of confidence 
in a student's t-distribution for rif — 23 — 2 = 21 
degrees of freedom. 

From the relationship between T a and -R max , 



T a = 79.7- 0.251i? n 



9.3, 



(5) 



one can estimate the rise time T a for cycle 24. 
Using the predicted value (84) of i? max (24), one 
obtains T a (24) = (59 ± 9) months. Therefore, the 
peak of cycle 24 may probably occur around Oc- 
tober, 2013, slightly later than that (May, 2013) 
by both NASA Marshall Space Flight Center 
(MSFCjl based on a quasi-Planck function and 
NOAA space weather prediction center (SWPCjl 
based on a consensus decision of "The Solar Cycle 
24 Prediction Panel, and that (J une, 2013) based 
on a modified Gaussian function (Du 2011b[) . 



4. Discussions and Conclusions 

Studying the correlation between -R max of a so- 
lar cycle and a related parameter is useful to un- 
derstand the dynamic process of the cycle. A high 
correlation can be use d to estimate the strength 
of a new solar c ycle (iKand 20101: IPesneH 120081: 
Messerotti et al.1 120091: |Pu| ^OllcL IWangetall 
20081 I2002L 120091 lLd l2004i riLil l2009l) . For exam- 
le, Ohl's geomagnetic precursor method (Ohl 



197G) succeeded in predicting -R max in cycles 20- 
22 (|Lavdenlll99lUThompson| [l993; Shastri 1998; 



2 http: / /solarscience. msfc.nasa.gov/SunspotCycle.shtml 
3 http: / /www. swpc.noaa.gov/SolarCycle/index. html 



Schussledl2007h due to the high correlation coef- 
ficients (> 0.8) between i? max and geomagnetic- 
based parameters. However, a high correlation 
does not always yield a satisfactory prediction 
(iDu et all l2009bt (Du fc Wand 1201 lat IPul 1201 la : 
Cameron fc Schiisslerl 120071: |Pu fc Wand 12008. 



2010h and a low correlation may also yield an 



curate prediction in some cases (see Section [272]) . 

A prominent feature in the solar cycle is the so- 
called Waldme ier effect that str o nger cycles tend 
to rise faster IWaldmeierl Il939t lHathaway et al 
2002; IDu et al.ll2009atlCameron fc Schussferll200 



This effect has already begun to work in the early 
phase of the cycle (/3 a ). The variation in /3 a reflects 
the temporal evolution of the dynamic process of 
solar magnetic activities from the initial phase to 
the peak phase of the cycle, and so, /3 a can provide 
very useful information for the cycle. 

In this study, we analyzed the temporal varia- 
tion in the correlation coefficient (r) between i? max 
and /3 a as a function of Am months after the so- 
lar minimum (i? m i n ) and the predictive power of 
(3a. on i? max as the solar cycle progresses. First, 
it is shown that r increases as Am increases with 
an inflection point over 0.8 at about A771 = 20 
months. The dynamic process of the solar activ- 
ity is more non-linear near the initial phase of the 
cycle (r < 0.5 if Am < 10) and tends to be stable 
after twenty months entering the cycle. 

Besides, (3 a behaved rather well in predicting 
i? max : the prediction error is usually within the 
estimated one after about Am = 20 months en- 
tering a solar cycle, \R p (Am, n) — i? max (n)| < 
t T (n — 1)ct(Ato, n — 1) at the 90% level of confi- 
dence. This is a crucial point in prediction be- 
cause a method will be less useful if the prediction 
is not within the prediction range derived from 
the method. In addition, the relative prediction 
error (E) based on /3 a is usually less than 20% 
when Am > 20 months. Thus, /3 a is a good indi- 
cator for the subsequent i? max . Finally, E tends 
to decrease as the cycle progresses. Therefore, the 
maximum amplitude of a new cycle (i? max ) can be 
well estimated at twenty months after the start. 

It should be noted in Fig. [2] that the correla- 
tion between i? max and /3 a is not strong near the 
initial phase of the cycle, while the prediction of 
Rma,x from /3 a is rather good in some cycles (15, 
20, 22 and 23). Therefore, a high correlation is not 
the sole condition to obtain a more accurate pre- 
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Table 2: Some predictions for cycle 24 

Author 



Kane (2010) 
Choudhuri ct al. 
NASA/MSFC 
Du ( 2011b') 
S valgaardetaL 
Schatten (2005 ) 
Li et al. (2005) 




current study 
Jiang et al. (2007 ) 
Li et al. (2011) 
NOAA/SWPC 
Wang fc Sheelev (2009^ 
Wang et al. (2009) 
Wang et al. (2002) 
Hirematn*72008) 
Rigozo et al. (2011) 
Dabas et al. (2008) 
Tlatov (2009) 
Du fc Wang (2008) 
Hathaway fc Wilson (2006) 
Dikpati fe QilmanT2006) 





x ± a 


58 


±25 


68 




70 




72 


±11 


75 


±8 


80 


±30 


80 


/137 


82 


/53 


84 


±17 


84 


±19 


85 




88 




90 




97 


±25 


100 


±8 


101 


±18 


110 


±11 


113 




124 


±23 


135 


±12 


140 


±16 


160 


±25 


165 


±15 



method or predictor 



aa minimum 

solar dynamo model 

quasi-Planck function 

modified Gaussian function 

polar field 

polar field 

slow/fast riser 

aa minimum/corrected 

similar cycles 

rising rate 

solar dynamo model 

sunspot minimum 

consensus 

open flux 

similar cycles 

similar cycles 

autorcgrcssion 

spectral components 

geomagnetic disturbed days 

G oc 1/R Z 

cycle length 

I component of aa 

flux-transport dynamo model 



dynamo mode ls have begun to be applied in pre- 
dicting i? max dChoudhuri et all 120071: Jiang et al 
120071 iDikpati fc Gilmanl^Ooi T the predictions of 



diction rtCameron fe Schiisslerl I2OO7I: IPu fc WansJ 
20081 l2010t ISvalgaard et al .1120051 : ISchattenll2005l h 

which may depend on the variat ion of the corre- 
lation or long-term periodicities dDu et al.ll20 09b; 
Du fc Wandl2011al iDullioilal IPu fc Wandl2008 . 
201dh . In this study,i? max can be well estimated 
from /3 a even if the correlation coefficient (r) is not 
strong near the initial phase in some cycles, e.g., 
E < 20% for small r (< 0.5) at small Am in cycles 
15, 20, 22 and 23 (see Fig.©. 

Based on /3 a at the current state Am — 27, the 
peak sunspot number of the ongoing cycle (24) 
is predicted to be i? max = 84 ± 33 at the 90% 
level of confidence or i? max = 84 ± 19 with a 1-a 
uncertainty. This prediction is higher than a few 
predictions and lower than many others, some of 
which are shown in Tabled] 

Accurately predicting the strength of an up- 
coming solar cycle is important for both solar 
physics and solar-terrestrial environment. A re- 
liable prediction of fl may may test mo dels for ex- 
plaining the solar cycle (jPesnelll 120081 ). So far, a 
great many results have been published on the pre- 
diction of -R m ax for cycle 24, of which some are 
based on statistics and some others are related to 
physics (see Table [2]) . As the solar activity near 
the minimum between cycles 23 and 24 lasts so 
long a time at a low level before rising ( as shown 



in th e most spotless days since cycle 16, iLi et al 
201l[ l2010h . cycle 24 is unusual, which is draw- 
ing greater attention than ever. Besides, as solar 



the strength of cycle 24 attract special attention 
in order to test the predictive skill of solar dynamo 
models. 

Discrepancies are found in the predictions of 
-Rmax for cycle 24 by erent methods (TableEJ. Our 
prediction (84) is near to those by the polar field 
(or solar dynamo model), about 30% lower than 
the peak size of cycle 23. Recently, we find that 
cycle 24 is most lik ely similar to cycles 14 and 10 
(|Du fc Wangll2011bl ). Therefore, even if cycle 24 is 
not a strong cycle, large eruption events may also 
occur as in cycle 10 for the largest solar storm of 
the year 1859 (Carrington Event). 

Conclusions are summarized below. 

1. The correlation coefficient (r) between the 
maximum amplitude (-R max ) of a solar cycle 
and the rising rate (/3 a ) at Am months after 
the solar minimum (i? m i n ) increases as Am 
increases with an inflection point at about 
twenty months entering the cycle. 

2. The prediction error based on the linear re- 
lationship between i? max and /3 a is usually 
within the estimated one when Am > 20, 
|i? p (Am, n)-R max (n)\ < t T (n-l)a(Am,n- 
1), where cr(Am, n— 1) is the standard devi- 
ation of the regression equation for the data 
up to cycle n — 1, and t T (n — 1) is the t- value 
at the 90% level of confidence in a student's 
t-distribution. 

3. The relative prediction error (E) from the 
above technique tends to decrease as the cy- 
cle progresses and will be less than 20% when 
Am > 20. 

4. The current cycle (24) is temporarily pre- 
dicted to peak around October 2013 with a 
size of i? max = 84 ± 33 at the 90% level of 
confidence. 
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